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Abstract 


We  report  on  -120  days  of  observations  of  Cygnus  X-1  with  OSSE 
onboard  the  Compton  Observatory.  Emission  is  detected  in  the  range  50 
keV  to  1  MeV  and  we  find  evidence  for  a  continuum  of  hard  X-ray  flux 
levels  rather  than  the  existence  of  distinct  flux  states.  Comparisons  of  the 
source  spectra  with  various  theoretical  models  show  that  an  exponentially 
truncated  power  law  best  describes  the  average  spectrum  in  the  OSSE 
energy  band.  Although  we  have  measured  a  new  minimum  in  the  hard  X- 
ray  flux  from  the  source,  no  evidence  was  found  for  either  a  broad  1  MeV 
feature  or  a  narrow  511  keV  line  previously  reported  in  association  with  a 
low  flux  state.  Upper  limits  on  such  emission  features  are  an  order  of 
magnitude  lower  than  earlier  reported  detections.  The  5.6-day  periodicity  of 
the  source  measured  at  optical  wavelengths  was  not  detected  with  a 
sensitivity  to  the  rms  modulation  fraction  of  5%  in  the  60-140  keV  energy 
band. 


Subject:  Gamma  Rays:  Observations-Black  Holes-stars:individual  (Cyg  X-1)- 

X-Rays:Binaries 


1  INTRODUCTION 


Cyg  X-1,  one  of  the  brightest  and  most  extensively  studied  gamma  ray  sources  in 
the  sky,  is  a  5.6-day  binary  system  consisting  of  a  blue  supergiant  (HDE  226868)  and  a 
compact  companion  of  mass  ~6.3  M  (Gies  et  al.  1986,  Dolan  et  al.  1989).  Since  the 
compact  object  has  a  mass  greater  than  3  M  ,  Cyg  X-1  is  generally  considered  to  be  a 
prime  black  hole  candidate,  and  it  has  become  the  prototype  for  back  hole  characterization 
and  modeling.  The  possible  presence  of  features  in  the  spectrum  of  Cyg  X-1  has  been 
suggested  as  a  useful  criterion  for  distinguishing  black  hole  systems  from  those  containing 
neutron  stars  (Liang  et  al.  1992).  In  particular,  features  such  as  soft  X-ray  excesses  or 
gamma-ray  bumps  in  the  hard  X-ray/gamma-ray  energy  band  have  been  proposed  as 
signatures  of  black  holes,  giving  special  importance  to  an  improved  knowledge  of  the 
detailed  shape  of  the  Cyg  X-1  spectrum  at  high  energies. 

The  high-energy  radiation  from  Cyg  X-1  is  thought  to  be  emission  from  an 
accretion  disk  surrounding  a  black  hole  (Shapiro,  Lightman  &  Eardly  1976).  The  soft  1- 
10  keV  X-ray  spectrum  has  traditionally  been  described  by  a  dominant  “low”  state  with 
occasional  periods  of  “high”  state  emission  (e.g.,  Priedhorsky,  Terrell  and  Holt  1983). 
The  hard  X-ray  flux  is  well  represented  by  a  thermal  Comptonization  spectrum  (Sunyaev 
and  Titarchuk  1980;  Ling  et  al.  1983).  Using  HEAO-3  data.  Ling  et  al.  (1987)  identified 
three  distinct  levels  for  the  45-140  keV  flux  occuring  during  the  low  X-ray  state,  which 
they  labeled  yl,  y2  and  y3,  with  yl  being  the  state  with  the  lowest  45-140  keV  flux.  The 
intensity  of  the  hard  X-ray  flux  may  be  anti-correlated  with  the  intensity  of  the  soft  X-ray 
flux  (Bassani  et  al.  1989)  .  A  gamma-ray  excess,  characterized  by  a  1.2  MeV  wide 
gaussian  at  1  MeV  with  a  line  flux  of  1.6  x  lO"^  photons  cm"^  s"l,  was  detected  in  fall 
1979  when  Cyg  X-1  was  in  the  yl  state  (Ling  et  al.  1987).  A  weak  2a  detection  of  a 
narrow  511  peak  was  later  reported  from  analysis  of  the  HEAO  data  during  the  same 
episode  (Ling  et  al.  1989).  Other  groups  have  also  reported  measurements  of  MeV 
emission  from  Cyg  X-1  (see  Owens  and  McConnell.  1992  for  a  review  of  the  MeV 
observations,  including  calculations  of  equivalent  line  flux  intensities).  Lor  example.  Baker 
et  al.  (1973)  described  a  spectral  feature  near  2  MeV,  with  positive  emission  up  to  6  MeV. 
Mandrou  et  al.  (1978)  observed  fluxes  up  to  3  MeV,  although  they  did  not  see  any  bumps. 
McConnell  et  al.  (1989)  reported  detecting  a  hard  spectral  component  out  to  9  MeV.  A 
general  review  of  Cyg  X-1  studies  was  written  by  Liang  and  Nolan  (1984). 


Although  there  is  a  wide  spread  of  energies  over  which  MeV  emissions  have  been 
reported,  these  measurements  suggest  periods  of  spectral  hardening  in  the  MeV  region 
beyond  a  simple  extrapolation  from  the  hard  X-ray  spectrum.  There  are,  however,  a 
considerable  number  of  experiments  (Schoenfelder  et  al.  1974;  White  et  al.  1980;  Nolan 
1982;  Varendorff  et  al.  1990)  reporting  featureless  hard  X-ray  and  gamma-ray  spectra,  and 
the  limits  derived  from  these  experiments  are  well  below  the  flux  values  measured  by  other 
experiments.  In  addition,  Harris  et  al.  (1993)  derived  upper  limits  for  flaring  episodes  of 
12  or  more  days  from  the  long-term  data  base  of  the  gamma-ray  spectrometer  on  SMM, 
and  McConnell  et  al.  (1994)  recently  reported  no  evidence  for  MeV  features  from  Cyg  X-1 
in  data  from  the  COMPTEL  experiment  on  the  Compton  Observatory. 

The  presence  of  enhanced  emission  near  1  MeV,  such  as  detected  with  HEAO-3, 
was  interpreted  by  Eiang  and  Dermer  (1988)  as  the  emergent  radiation  from  a  photon- 
starved  pair-dominated  thermal  plasma.  In  this  picture,  the  hard  X  rays  are  produced  by 
thermal  Comptonization  of  a  copious  soft  photon  source  farther  out  in  the  disk,  in  accord 
with  earlier  models  of  the  hard  X-ray  spectrum  of  Cyg  X-1  and  the  innermost  region, 
shielded  from  the  soft  photon  flux,  reached  MeV  temperatures.  Spectral  complexity  is  also 
found  in  the  X-ray  regime.  Done  et  al.  (1992)  demonstrated  the  presence  of  a  reflection 
component  in  EXOSAT  and  HEAO-1  observations,  and  argued  that  a  multi-temperature 
thermal  Comptonization  model  was  required  to  give  a  good  fit  to  the  Cyg  X- 1  spectrum. 
Recently,  Haardt  et  al.  (1993)  suggested  that  thermal  Comptonization  in  an  accretion-disk 
corona  is  simply  described  by  a  power  law  with  an  exponential  cutoff.  Zdziarski  et  al. 
(1994)  combined  optically  thin  Comptonization  and  reflection  into  a  model  consisting  of  an 
exponentially  cutoff  power  law  with  reflection  in  a  form  useful  for  data  analysis.  The 
excellent  sensitivity  of  the  OSSE  measurements  provide  spectra  with  adequate  precision  to 
test  these  theoretical  models. 

In  this  paper,  we  present  results  of  measurements  of  Cygnus  X-1  using  the 
Oriented  Scintillation  Spectrometer  Experiment  (OSSE)  on  the  Compton  Gamma  Ray 
Observatory  (CGRO).  Given  the  importance  of  testing  the  models  of  the  source  and 
resolving  the  conflicting  nature  of  some  of  the  previous  measurements,  OSSE  devoted  122 
days  to  observing  the  source  at  different  intensity  levels  during  the  first  four  years  of  the 
CGRO  mission.  The  shape  of  the  spectrum,  the  search  for  narrow  and  broad  spectral 
features,  the  light  curve  of  daily  flux  averages,  the  flux  levels  of  the  source,  and  correlations 
of  hardness  ratios  are  addressed.  We  describe  OSSE  observations  of  Cyg  X-1  in  Section 
2,  the  temporal  variability  analysis  in  Section  3,  and  spectral  analysis  in  Section  4.  Our 


search  for  spectral  features  is  presented  in  Section  5,  and  discussion  of  the  results  is  given 
in  Section  6. 

2  OBSERVATIONS 

The  Oriented  Scintillation  Spectrometer  Experiment  is  one  of  four  instruments  on 
the  Compton  Gamma  Ray  Observatory.  The  observatory  was  launched  on  board  the  space 
shuttle  Atlantis  on  5  April  1991  and  has  been  operating  in  a  28  degree  inclination  low  earth 
orbit.  OSSE  consists  of  four  identical  and  independently  moveable  NaTCsI  phoswich 
detectors  sensitive  to  gamma  rays  between  50  keV  and  10  MeV.  Each  detector  has  a 
tungsten  collimator  which  defines  a  3.8°  x  11.4°  EWHM  rectangular  field  of  view,  and 
active  shielding  provided  by  Nal  anticoincidence  detectors  configured  in  an  annulus  around 
the  main  detector.  Each  detector  has  an  effective  area  of  ~  500  cm^  at  5 1 1  keV.  Eight-  or 
sixteen-second  spectra  are  accumulated  on  board  and  then  transmitted  to  Earth  for 
spectroscopic  analysis,  while  limited  event-by-event  or  fast  rate  samples  over  certain 
energy  bands  are  also  accumulated  in  parallel  for  timing  analysis.  We  will  report  here  on 
analysis  of  the  spectroscopy  data.  Timing  analysis  of  these  data  at  mHz  to  kHz  frequencies 
is  not  considered  here  and  will  be  presented  elsewhere  (see  Bridgman  et  al.  1994  for  an 
early  analysis).  Spectra  are  integrated  into  two  minute  intervals,  alternately  viewing  the 
source  field  and  background  fields  offset  by  ~5  degrees.  Using  the  known  response  of  the 
detectors,  the  background-subtracted  count  spectra  are  converted  into  photon  spectra  using 
a  forward-folding  technique.  A  full  description  of  the  instrument  and  standard  methods  of 
data  analysis  are  given  by  Johnson  et  al.  (1993). 

OSSE  observed  Cyg  X-1  during  17  intervals,  each  spanning  one  or  more  days  in 
the  first  four  years  of  the  CGRO  mission,  for  a  total  of  ~2.4  x  10^  live  seconds  of 
background-subtracted  spectroscopy  data.  Table  1  summarizes  these  observations.  The 
large  spread  in  total  live  times  of  the  different  observations  resulted  from  variations  in  the 
number  of  detectors  used  in  the  observation  (1  to  4)  and  the  length  and  quality  of  each 
observation.  These  data  are  not  a  random  sample  of  the  flux  levels  of  the  source  because 
four  of  the  observations,  numbers  6,  7,  8  and  10,  were  in  response  to  a  low  hard  X-ray 
state  detected  by  the  Burst  and  Transient  Source  Experiment  (BATSE)  on  CGRO 
(Eishman  et  al.  1989).  The  observations  were  analyzed  for  each  detector  and  each  day 
separately,  and  the  spectra  were  added  at  the  end  of  the  analysis. 


3  TEMPORAL  ANALYSIS 


The  integrated  flux  in  the  45-140  keV  band  for  each  day  of  data  is  shown  in  Figure 
1,  where  the  x-axis  represents  days  of  data  with  a  short  break  between  viewing  periods  that 
is  not  representative  of  the  actual  time  between  the  observations.  The  energy  band  was 
chosen  to  match  the  band  used  in  the  analysis  of  the  HEAO-3  data.  The  data  between  60 
and  140  keV  were  fit  with  a  thermal  bremsstrahlung  spectrum  and  the  resulting  fit  was 
used  to  calculate  a  45-140  keV  flux.  For  the  days  where  the  instrument  configuration 
provided  spectra  down  to  40  keV,  the  calculated  and  measured  flux  were  in  agreement. 

As  seen  in  Figure  1,  there  are  no  obvious  discrete  states,  but  rather  a  continuum  of 
intensities.  This  becomes  clearer  when  the  daily  flux  is  observed  over  a  long  period  of 
time,  as  is  done  by  the  BATSE  experiment  (Paciesas  et  al.  1994).  Compared  with  the 
historic  yl,  y2  and  y3  levels,  the  observations  of  Eebruary  1994  (CGRO  viewing  period 
318)  revealed  a  flux  intensity  of  ~  1/5  the  previously  reported  lowest  flux  in  this  band.  The 
source  has  therefore  a  greater  range  of  intensities  than  previously  measured.  Note, 
however,  that  no  OSSE  measurements  have  been  performed  at  intensities  comparable  to 
the  y3  level. 

The  light  curve  shown  in  Eigure  1  can  be  searched  for  periodicities,  in  particular  for 
the  5.6-day  orbital  periodicity.  The  data  set  is  sparse,  with  only  -100  days  of  data  out  of  a 
1024  day  interval,  with  the  contiguous  daily  data  sets  varying  from  two  days  to  two  weeks 
in  length.  The  measurements  are  therefore  not  very  sensitive  to  periodicities  between  a  few 
weeks  and  a  few  months.  Eor  this  analysis  the  data  were  binned  into  individual  days; 
therefore  the  shortest  detectable  variability  is  two  days.  We  calculated  the  East  Eourier 
Transform  of  the  mean  subtracted  data  set  for  the  first  1024  days  of  the  mission  and  found 
most  power  at  low  frequencies  due  to  overall  trends  in  the  data;  more  specifically,  the 
significant  drop  in  flux  in  Eebruary  1994.  To  detrend  the  data,  we  subtracted  a  fifth  order 
polynomial  in  time  and  calculated  the  East  Eourier  Transform  of  the  residuals.  We  find  no 
evidence  for  periodicity  in  this  sparse  data  set,  other  than  a  peak  at  very  low  frequency  due 
to  the  data  sampling.  The  upper  limit  to  the  amplitude  of  a  5.6-day  periodicity  is  5xl0‘3 
photons  cm‘2  s‘l  in  the  60-140  keV  band,  corresponding  to  an  rms  modulation  fraction  of 
less  than  5%. 


4  SPECTRAL  ANALYSIS 


The  spectrum  for  the  sum  of  all  data  through  May  1995  is  showu  iu  Figure  2.  The 
spectrum  is  smooth,  without  breaks  or  liue  features.  It  exteuds  from  60  keV,  the  lower 
threshold  used  iu  this  spectral  aualysis,  to  ~1  MeV,  where  the  flux  level  falls  below  the 
OSSE  seusitivity.  The  high  statistical  precisiou  of  the  Cyg  X-1  spectrum  permits 
differeutiatiou  amoug  theoretical  models  to  much  higher  euergies  thau  previously  possible. 
We  use  the  four  differeut  fuuctioual  forms  to  fit  the  data.  The  first  is  a  simple  two- 
parameter  thermal  bremsstrahluug  model  described  by 

d>(E)  =  (1) 

where  d>(E)  is  the  flux  at  euergy  E,  K  is  the  uormaliziug  coustaut  aud  T  is  the  temperature. 
The  secoud  fuuctioual  form  we  use  is  the  three-parameter  thermal  Comptouizatiou 
spectrum  of  Suuyaev  aud  Titarchuk  (1980).  It  is  derived  for  a  spherical,  uuiform  electrou 
plasma  with  Thomsou  depth  x  =  UeOxRj  where  Ug  is  the  electrou  deusity  aud  R  is  the 
radius  of  the  regiou,  aud  is  strictly  valid  iu  the  uourelativistic  (kT  «  mgC^),  diffusive  (x 
»  1)  regime.  Higher  temperature  correctious  have  beeu  derived  by  Titarchuk  (1994),  aud 
Hua  aud  Titarchuk(1995),  but  Moute  Carlo  simulatious  iucludiug  pair  productiou  (Skibo  et 
al.  1995)  are  required  to  provide  accurate  represeutatious  at  all  temperatures  aud  optical 
depths.  However,  the  Moute  Carlo  results  aud  the  refiued  model  of  Titarchuck  (1994) 
agree  closely  with  the  origiual  model  of  Suuyaev  aud  Titarchuk  (1980)  for  the  parameter 
values  of  Cyg  X- 1 . 

The  third  fuuctioual  form  employed  is  a  three-parameter  expoueutially  truucated  power  law 
giveu  by 


=  (2) 

where  E  is  the  photou  spectral  iudex  aud  Eg  is  the  cutoff  euergy.  We  also  iuclude  the 
possibility  that  some  fractiou  of  the  emissiou  described  by  equatiou  (2)  is  iucideut  ou  a 
cold  optically  thick  scatteriug  medium  for  the  fourth  fuuctioual  form  we  employ.  The 
effects  of  Comptou  reflectiou  are  coutaiued  iu  a  code  provided  by  Zdziarski  (1995a).  The 
model  coutaius  a  reflectiou  parameter  that  is  proportioual  to  the  solid  augle  covered  by  the 
reflectiou  medium:  a  reflectiou  coefficieut  of  1  is  a  source  just  above  au  iufiuite  plaue  of 
reflectiug  material,  a  coefficieut  of  two  implies  that  the  sources  is  embedded  iu  the 
reflectiug  medium. 


Table  2  lists  the  fits  of  these  various  models  to  the  sum  of  all  the  Cyg  X-1  data. 
The  normalizing  energy  was  fixed  at  100  keV  and  all  daily  fits  were  over  the  energy  range 
from  60  keV  to  4  MeV.  When  all  the  observations  are  added  together,  the  systematic 
differences  among  detectors  become  important  and  the  spectra  have  to  be  fit  detector  by 
detector  with  the  various  models.  The  difference  between  the  variance  in  the  distribution  of 
the  parameters  obtained  and  the  square  of  the  error  in  the  fit  parameters  is  used  as  a 
measure  of  systematic  errors.  The  statistical  errors  are  68%  confidence  intervals.  The 
parameters  in  Table  2  are  listed  as  the  average  of  that  parameter  for  all  detectors  with  the  fit 
error  followed  by  the  systematic  error,  in  parentheses.  The  reduced  quoted  is  the 
average  per  degree  of  freedom  of  the  fits  over  the  detectors.  The  best  fit  models  given 
by  the  functional  form  of  equation  (2)  and  a  thermal  Comptonization  model  are  overlaid  on 
the  data  in  Figure  2. 

The  model  with  an  exponentially  truncated  power  law  gives  the  best  fit  to  the  data. 
The  thermal  bremsstrahlung  and  the  thermal  Comptonization  models  both  diverge  from 
the  measured  spectrum  at  higher  energies  where  they  underestimate  the  flux  from  the 
source.  The  exponentially  truncated  power  law  with  reflection  produces  a  fit  to  the  data 
slightly  worse  than  without  reflection,  and  the  fit  improves  as  the  amount  of  reflection 
decreases.  The  difference  in  the  photon  spectra  between  the  model  with  and  without 
reflection  occurs  mostly  at  the  lower  end  of  the  OSSE  energy  range,  where  the  systematic 
errors  in  the  instrument  efficiency  and  calibration  are  largest.  The  confidence  in  rejecting 
models  based  solely  on  the  lowest  energy  bins  is  therefore  smaller  than  the  confidence  in 
rejecting  models  which  diverge  at  higher  energies  like  the  thermal  bremsstrahlung  and 
thermal  Comptonization  models.  The  exponentially  truncated  power  law  is  also  the  best 
model  when  each  observation  is  summed  and  fitted  separately.  Table  3  lists  the  results  of 
the  fits  to  each  observation,  along  with  the  source  luminosity  in  the  60  keV-1  MeV  range  , 
assuming  a  Cyg  X-1  distance  of  2.5  kpc.  The  statistical  error  in  the  luminosity  is  1%  or 
better  for  all  observations. 

The  optically  thin  thermal  bremsstrahlung  model  still  provides  a  good  fit  to  the 
daily  spectra.  This  simple  two-parameter  model  is  used  to  study  the  temperature-intensity 
dependence  of  the  source.  Figure  3  shows  that  the  temperature  of  the  source  and  the 
intensity  vary  within  a  limited  range,  except  for  a  few  excursions  to  the  low  temperature- 
low  amplitude  regime.  During  the  very  low  flux  observation  of  February  1994,  the 
temperature  of  the  source  was  ~  70  keV  compared  with  110  keV  for  the  other 
observations.  Note  that  the  scatter  in  the  points  on  this  figure  and  the  following  figures  is 


much  larger  than  the  statistical  or  systematic  errors.  Alternatively  fitting  the  daily  data  with 
the  thermal  Comptonization  model  derived  by  Sunyaev  and  Titarchuck  (1980),  we  find  a 
similar  scatter  in  the  optical  depth  parameter,  with  the  temperature  fairly  constant  with  a 
mean  of  65  keV  and  a  standard  deviation  of  less  than  5  keV.  With  the  temperature  fixed  at 
65  keV,  the  optical  depth-intensity  distribution  is  as  shown  in  Figure  4. 

In  Figure  5,  we  display  spectra  in  terms  of  E2d>(E)  in  order  to  highlight  the  energy 
of  peak  power  output  in  the  OSSE  band  and  emphasize  the  differences  between  spectra 
measured  at  different  flux  states.  The  spectrum  with  open  squares  is  the  sum  of  all  days 
with  a  flux  (in  the  45-140  keV  band)  of  less  than  0.06  photons  cm‘2  s'b  The  spectra  with 
filled  squares,  open  circles  and  filled  circles  are  for  a  flux  between  0.06  and  0.08,  between 
0.08  and  0.105  and  greater  than  0.105  photons  cm-2  s‘l,  respectively.  The  low  amplitude 
spectra  are  more  like  a  power  law  than  the  other  data.  If  a  power  law  is  assumed  for  the 
observation  with  the  lowest  amplitude  spectrum,  observation  10,  the  best  fit  is  for  an 
amplitude  of  7.5  ±  0.1  x  10‘2  photons  cm‘2  s‘l  MeV‘l  at  100  keV  and  a  power  law  index 
of  -2.72±.03  with  a  probability  of  0.3. 

5  SEARCH  FOR  LINE  FEATURES 

To  search  for  any  spectral  features  associated  with  specific  flux  levels  that  would 
have  been  averaged  out  in  the  total  spectrum,  the  data  were  divided  into  individual  days  and 
each  day  was  analyzed  separately.  Following  the  analysis  of  Eing  et  al.  (1987),  we 
searched  for  the  presence  of  MeV  bumps  by  fitting  each  day  of  data  with  a  thermal 
bremsstrahlung  spectrum  and  a  1.2  MeV  wide  gaussian  line  feature  centered  at  1  MeV. 
Near  the  yl  flux  level,  an  increase  in  the  fit  1  MeV  amplitude  would  be  expected.  The  fit 
amplitudes  are  shown  in  Figure  6a.  We  see  no  significant  amplitude  for  a  broad  gaussian 
feature  for  any  flux  level  of  the  source.  When  all  the  data  are  added  together,  as  in  Figure 
2,  the  90%  confidence  level  (CE)  upper  limit  to  the  1  MeV  line  is  4  x  10""^  photons  cm‘2  s" 
1,  compared  to  the  flux  of  1.6  x  10"2  photons  cm‘2  s‘l  measured  by  Eing  et  al.  The  95% 
CE  upper  limit  for  a  1-1.5  MeV  wide  feature  near  4  MeV  is  2.0  x  lO"^  photons  cm‘2  s‘l  . 

The  search  for  lines  can  be  repeated  at  any  energy  and  width.  Given  the  previous 
reports  of  a  narrow  511  keV  line,  we  fit  the  data  for  a  narrow  (<  instrument  resolution)  511 
keV  line  for  each  day  of  data.  These  results  are  shown  in  Figure  6b.  We  do  not  find  any 
days  with  a  significant  excess  511  keV  emission  nor  a  correlation  between  the  511  keV 
emission  and  the  overall  source  strength.  The  90%  CE  upper  limits  to  the  amplitude  are 


between  0.5  and  1.0  x  10  ^  photons  cm‘2  s‘l  for  most  days.  The  90%  CL  upper  limit  for  a 
narrow  511  keV  in  the  total  spectrum  is  7  x  10"^  photons  cm‘2  s‘l. 

Since  Ramaty  et  al.  (1981)  predict  a  broad  511  keV  feature  (FWHM  -160  keV), 
we  also  searched  for  a  160  keV  wide  peak  at  511  keV.  The  results  show  no  significant 
positive  detections.  The  fit  flux  values  for  the  broad  line  are  significantly  larger  than  for  the 
narrow  line  because  the  broad  line  has  a  larger  effect  on  the  underlying  continuum 
spectrum  represented  by  a  thermal  bremsstrahlung  model.  Still,  the  derived  flux  values 
provide  new  limits  to  theoretical  models  with  a  90%  CL  upper  limit  of  1  to  3  x  10"^ 
photons  cm‘2  s‘l  for  most  days.  The  90%  CL  upper  limit  for  a  broad  511  keV  feature  in 
the  total  spectrum  is  2  x  10""^  photons  cm‘2  s‘C 

6  DISCUSSION 

The  variability  of  Cyg  X-1  at  gamma  ray  energies,  as  seen  Figure  1,  is  a  factor  of 
-15,  much  greater  even  than  the  factor  of  -3  measured  during  the  HEAO-3  mission. 
However,  the  variability  does  not  seem  to  be  between  discrete  states  but  rather  among  a 
continuous  range  of  possible  flux  values.  The  optically  thin  thermal  bremsstrahlung  fits  to 
the  individual  days  of  data  show  that  the  source  usually  fluctuates  in  a  seemingly  random 
fashion  in  a  narrow  range  of  intensity-temperature  phase-space.  This  is  unlike  the 
behaviour  observed  in  GRO  J0422-I-32,  where  the  temperature  is  monotonically  increasing 
in  time  as  the  transient  decreases  in  brightness  (Kroeger  et  ah,  1995). 

The  spectrum  is  best  described  by  a  model  with  optically  thin  Comptonization. 
However,  a  preliminary  analysis  reveals  that  this  model  cannot  be  applied  to 
contemporaneous  GINGA  (Ebisawa,  1994)  and  OSSE  data.  A  more  complex  model 
using  an  exponentially  truncated  power  law  with  an  angle  dependent  reflection  component 
might  fit  both  sets  of  data  simultaneously  (Zdziarski  et  ah,  1995b).  Using  the 
exponentially  truncated  power  law  model  to  find  the  maximum  in  the  luminosity  per 
decade  and  correlating  that  with  the  luminosity  in  the  60  keV-1  MeV  band,  the  relationship 
shown  in  Eigure  7  is  obtained.  The  maximum  in  the  distribution  of  luminosity  per  decade 
seems  correlated  to  the  total  luminosity.  The  spectrum  for  the  lowest  intensity  observation, 
observation  10,  is  then  interpreted  as  the  same  model  with  a  maximum  at  «  60  keV  in  the 
distribution  of  luminosity  per  decade,  instead  of  the  more  typical  100  keV,  causing  the 
power  law-like  appearance  of  this  spectrum. 


Extrapolating  the  exponentially  truncated  power-law  model  to  lower  energies,  we 
find  that  the  calculated  flux  in  the  3-11  keV  range  for  observation  10  (the  lowest  intensity 
observation  in  the  OSSE  energy  band)  is  ~3  times  larger  than  for  observation  1,  a  typical 
observation.  This  anticorrelation  between  the  extrapolated  low  energy  X-ray  flux  and  the 
measured  gamma  ray  flux  is  consistent  with  earlier  obervations  (Priedhorsky  et  al.  1983) 
and  is  not  very  model  dependent.  The  crossover  energy  for  the  low  intensity  gamma  ray 
spectra  to  become  dominant  is  between  10-20  keV,  depending  on  which  spectra  are 
compared. 

We  have  searched  for  the  amplitude  of  a  broad  MeV  bump,  a  narrow  511  keV  line, 
and  of  a  broad  511  keV  line.  The  distribution  of  fit  MeV  amplitudes  (Eigure  8)  is 
gaussian,  with  a  small  non-zero  mean,  indicating  the  underlying  thermal  bremsstrahlung 
distribution  under-estimates  the  continuum  flux  at  1  MeV  by  ~5  x  10""^  photons  cm‘2  s‘l. 
There  is  no  evidence  for  an  excess  of  points  at  high  flux  values,  and  no  single  day  with  an 
amplitude  near  that  measured  by  Eing  et  al.  (1987)  over  a  12  day  period.  The  width  of  the 
distribution  in  Eigure  8  is  consistent  with  the  uncertainty  in  the  amplitudes  histogrammed. 
Since  no  strong  signal  was  found,  upper  limits  are  derived  for  such  features.  These  upper 
limits  are  lower  than  previous  measurements  and  constrain  the  models  for  Cyg  X-1.  The 
90%  confidence  level  upper  limit  to  the  1  MeV  bump  flux  is  in  the  1-4  x  lO"^  photons  cm" 
2  s"l  range  for  most  days.  This  clearly  demonstrates  that  the  MeV  bump  is  not 
systematically  associated  with  the  yl  state.  If  MeV  bumps  of  the  intensity  reported  by 
Eing  et  al.  do  exist,  they  must  be  transient  within  the  yl  state.  Our  non-detection  of  MeV 
features  limits  the  duty  cycle  of  such  a  process.  Eor  1-day  emission  episodes,  there  is  a 
90%  probability  we  would  have  detected  a  duty  cycle  greater  than  2.3%,  or  more  than  25 
outbursts  over  the  first  three  years  of  the  mission.  Eor  4-day  episodes  the  duty  cycle  must 
be  less  than  4.5%,  and  for  14-day  episodes,  the  duty  cycle  must  be  less  than  9%,  which 
corresponds  to  ~6.5  outbursts  over  the  mission  so  far.  Harris  et  al.  (1994)  obtained  a  95% 
upper  limit  of  4.1%  on  the  duty  cycle  for  12-day  episodes. 

According  to  the  models  of  Dermer  et  al.  (1989)  and  the  compilation  of  data  by 
Bassani  et  al.  (1989),  there  should  be  an  anticorrelation  between  the  hard  X-ray  flux  and  the 
gamma  ray  flux.  We  have  binned  our  data  in  two  bands,  0.06-0.4  MeV  and  0.4- 1.5  MeV 
and  plotted  the  high  energy  band  versus  the  low  energy  band,  for  a  direct  comparison  with 
Bassani  et  al.  (Eigure  9).  Erom  the  total  data  set,  one  would  conclude  that  there  is  a 
correlation,  not  an  anticorrelation,  between  the  gamma  ray  flux  and  the  hard  X-ray  flux. 
This  is  consistent  with  our  earlier  interpretation  of  higher  temperatures  at  higher  flux.  The 


anticorrelation  modeled  by  Dermer  and  Liang  was  for  a  fixed  accretion  rate.  However,  it  is 
likely  that  the  changes  measured  are  due  to  a  change  in  accretion  rate.  In  that  case,  the 
predicted  hard  X-ray  and  gamma  ray  components  would  rise  and  fall  together  and  would 
be  correlated,  as  observed.  The  correlation  coefficient  for  these  data,  when  weighed  by  and 
including  only  the  error  in  the  high  energy  flux,  is  0.72,  corresponding  to  a  probability  that 
the  variables  are  not  correlated  of  <  0.01.  If  only  data  with  hard  X-ray  flux  >  0.20  photons 
cm‘2  s‘l  are  considered,  the  correlation  coefficient  of  this  ratio  is  0.03,  implying  no 
correlation  between  the  fluxes  during  most  observations. 

The  0.4  -1.5  MeV  band  is  not  optimized  for  daily  OSSE  observations  of  Cyg  X-1, 
since  the  source  is  not  detected  beyond  -800  keV  on  a  daily  basis.  This  energy  band  was 
selected  to  provide  a  direct  comparison  with  the  analysis  of  Bassani  et  al.  To  better  check 
our  compatibility  with  the  model  of  Dermer  et  al.  that  postulates  a  pivot  point  near  100 
keV,  we  have  derived  the  hardness  ratio  between  the  200-400  keV  band  and  the  60-80  keV 
band.  A  spectral  shape  independent  of  amplitude  would  produce  a  horizontal  line  in  this 
plot,  and  the  expected  anticorrelation  would  be  seen  as  a  decreasing  curve.  From  Figure 
10,  we  see  the  same  pattern  as  in  the  temperature- amplitude  scatter  plots.  There  is  a 
positive  correlation  (correlation  coefficient  of  0.32)  between  the  hardness  ratio  and  the  60- 
80  keV  flux  for  the  low  amplitude  points  only  (<  0.7  photons  cm‘2  s‘l  in  the  60-80  keV 
band),  again  consistent  with  higher  temperatures  at  higher  flux.  But  for  the  dense,  high 
hardness  ratio  cluster  of  data  (>  0.7  photons  cm‘2  s‘l  in  the  60-80  keV  band  and  >  0.035 
photons  cm‘2  s‘l  in  the  hardness  ratio),  there  is  an  anticorrelation  (correlation  coefficient  of 
-0.48)  between  the  ratio  and  the  low  energy  flux. 

We  hope  to  shed  further  light  on  this  correlation-anticorrelation  pattern  by  acquiring 
more  data  on  Cyg  X-1  when  it  is  at  intermediate  or  very  bright  flux  levels. 
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Table  1:  OSSE  observations  of  Cyg  X-1. 


No. 

Dates(^) 

Viewing 

Period 

Total  Live 

Seconds^^^)  xlO"^ 

1 

91/150-165 

2 

39.2 

2 

91/220-227 

7.0 

10.9 

3 

91/332-346 

15 

25.4 

4 

92/336-357 

203 

21.8 

5 

93/040-053 

209 

20.1 

6 

93/080-082 

212 

5.4 

7 

93/133-134 

221 

2.5 

8 

93/151-154 

223 

7.0 

9 

93/274-277 

303.4 

7.9 

10 

94/032-039 

318.1 

43.8 

11 

94/144-145 

328 

2.7 

12 

94/158-159 

331 

2.8 

13 

94/165-169 

331.5 

12.1 

14 

94/174-175 

332 

1.3 

15 

94/192-193 

333 

2.5 

16 

95/130-143 

419.5 

9.7 

17 

95/143-156 

420 

8.8 

(a)  dates  are  year/day  of  year. 

(b)  GRO  mission  timeline  viewing  period  identification. 


(c)  the  total  livetime  is  the  sum  of  the  livetimes  of  the  detectors  observing  the  source. 


Table  2:  Fits  of  various  models  to  the  122-day  summed  Cyg  X-1  spectrum. 


Model 

Amp.  @  100  keV 
(y  cm‘2  s‘l  MeV‘l) 

kT,  Ec 
(keV) 

Photon  Index  x 

%^/DOF 

Thermal 

0.476±0.002(±0.017) 

108±1(±5) 

50.0 

Bremsstrahlung 

Thermal 

0.465±0.001  (±0.021) 

65±1(±3) 

1.94±.04(±.05) 

19.6 

Comptonization 

Power  law  times 

0.470±0.001(±0.017) 

158±3(±12) 

1.39±0.02(±0.08) 

2.9 

exponential 

Power  law  times 

0.473±0.001(±0.017) 

156±3(±4) 

0.96±0.04(±0.08) 

8.3 

exponential 
with  refl.  @1.5 


(a)  For  the  Thermal  Bremsstrahlung  and  Sunyaev-Titarchuck  model,  the  parameter  is  kT. 

For  the  exponentially  truncated  power  law  models,  the  parameter  is  the  cutoff  energy,  Eq. 


Table  3:  Fits  of  an  exponentially  truncated  power  law 
to  the  individual  Cyg  X-1  observations. 


Observation 

Number 

Amp.  @  100  keV 
(y  cm‘2  s‘l  MeV‘l) 

X  10-3 

Ec 

(keV) 

Photon  Index 

Luminosity^) 

(erg  s-1) 

X  1037 

1 

624±1 

148±2 

-I.22±.0I 

1.30 

2 

657±3 

151±9 

-I.2I±.05 

1.39 

3 

826±2 

139±3 

-I.I8±.02 

1.68 

4 

752±1 

137±2 

-I.28±.02 

1.46 

5 

644±1 

150±3 

-I.34±.02 

1.28 

6 

777±7 

177±25 

-I.55±.0I 

1.54 

7 

613±3 

135±6 

-I.I4±.04 

1.24 

8 

528±2 

153±6 

-I.29±.03 

1.08 

9 

459±2 

155±7 

-I.70±.04 

0.81 

10 

76±1 

582±250 

-2.52±.08 

0.15 

11 

756±3 

138±6 

-I.30±.04 

1.46 

12 

499±3 

182±24 

-I.79±.09 

0.92 

13 

236±4 

299±67 

-2.I6±.09 

0.45 

14 

497±4 

188±25 

-I.64±.09 

0.98 

15 

516±4 

129±13 

-I.43±.09 

0.92 

16 

698±3 

137±8 

-I.I6±.05 

1.42 

17 

715±4 

124±7 

-I.0I±.06 

1.45 

(a)  Integral  luminosity  in  the  60  keV-lMeV  band. 
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Figure  1:  Flux  in  the  45-140  keV  band  for  each  day  of  Cyg  X-1  data.  The  line 
segments  show  the  viewing  period  associated  with  the  data  points. 

Also  shown  are  the  historic  yl  and  y2  states  as  defined  by  Ling  et  al.  (1987). 
Alternate  viewing  periods  are  shown  with  open  and  filled  squares. 
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Figure  2:  Spectrum  of  Cyg  X-1  for  all  the  data.  The  solid  line  is  the  best  fit  to  the 
data  with  an  exponentially  truncated  power  law  model  while  the  dotted  line  is  the  best  fit 
using  the  thermal  Comptonization  model  of  Sunyaev  and  Titarchuck  (1980).  Upper  limits 
are  2  sigma.  See  Table  2  for  the  parameters  of  the  fits. 
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Figure  5:  Photon  spectrum  multiplied  by  the  energy  squared.  The  data  were 
grouped  into  four  bins  according  to  their  45-140-keV  flux:  <0.06  photons  cm‘2  s‘l, 
between  0.06  and  0.08  photons  cm‘2  s‘l,  between  0.08  and  0.105  photons  cm‘2  s‘l,  and 
greater  than  0.105  photons  cm‘2  s‘l. 
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Figure  6:  Amplitude  and  error  of  the  1  MeV  gaussian  fit  (Figure  6a)  and  the 

narrow  511  keV  line  (Figure  6b)  on  top  of  a  thermal  bremsstrahlung 
continuum  for  each  day  of  data.  The  X-axis  is  the  same  discontimuous 
time  axis  as  in  figure  1 . 
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Figure  7:  Peak  in  luminosity  per  decade  versus  luminosity  between  60  and  1000 
keV  for  the  model  fits  of  all  the  observations,  as  in  Table  3. 
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Figure  8:  Distribution  of  the  daily  amplitudes  of  the  1  MeV  gaussians.  The 
amplitude  reported  by  Ling  et  al.  (1987)  is  shown  for  comparison. 


Gamma  Ray  Flux  (cm 


'c/D 


cq 


0.006 

0.004 

0.002 

0.000 


-0.002 


-0.004 


0.00  0.10  0.20  0.30 

Hard  X-ray  Flux  (cm'^  s'^) 


Figure  9:  Flux  in  the  0.4- 1.5  MeV  band  versus  the  flux  in  the  0.06-0.4  MeV 

band  for  each  day  of  data.  These  bands  are  nearly  identical  to  those  in 
the  compilation  of  Bassani  et  al.  (1989). 
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Figure  10:  Ratio  of  the  flux  in  the  200-400  keV  band  over  the  flux  in  the  60-80 
keV  band  versus  the  flux  in  the  60-80  keV  band  for  each  day  of  data. 


